TVion dynamics in coupled double quantum wells. Electron density effects. 
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We have studied the coherent dynamics of injected electrons when they are either free or bounded 
both in excitons and in trions (charged excitons). We have considered a remotely doped asymmetric 
double quantum well where an excess of free electrons and the direct created excitons generate 
trions. We have used the matrix density formalism to analyze the electron dynamics for different 
concentration of the three species. Calculations show a significant modification of the free electron 
inter-sublevel oscillations cWe have studied the coherent dynamics of injected electrons when they 
are aused by electrons bound in excitons and trions. Based on the present calculations we propose 
a method to detect trions through the emitted electromagnetic radiation or the current density. 

PACS numbers: 73.40.Gk, 78.47,+p 
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I. INTRODUCTION 

Coherent oscillations (quantum beats) of the electronic charge in double quantum wells (DQW's) have been widely 
studied both experimentally and theoreticalljiiiSi&iSiSiLS. The most usual method to achieve charge generation in these 
structures is the photoexcitation by an ultrashort laser pulse, which causes the same concentration of photoexcited 
electrons and holes and so, the corresponding excitons. Another possibility consists on injecting electrons and holes 
by doping two regions close to the structur e) 9 ! 10 . Injection offers the advantage of preventing the interaction between 
the electromagnetic field associated to the excitation and the excited electrons. Here excitons are directly created by 
hole-assisted electron resonant tunneling (see Fig. 1). Another characteristic of this type of excitons (the so-called 
direct-created excitons) is that they have spatial coherence, with an in-plane momentum k ea , c ii ~ 0, and can radiate 
in the direction perpendicular to the DQW. Thus, direct-created excitons are good candidates for their application in 
vertical planar microcavities. Moreover, injection has two additional advantages: the control of the electron and hole 
densities by means of the n (donor) or p (acceptor) impurity concentration plus an external applied electric fieldii, 
and the avoidance of exciton thermalization. This thermalization is an inevitable consequence of non-tuned laser 
excitation. 

Recently, trions (charged excitons) have been observed in photoexcited doped semiconductor quantum wells. The 
doping is necessary to get an excess of electrons. These electrons together with the excitons lead to the generation 
of the negatively charged excitonsiiASii^. These bound complexes of three particles have a binding energy large 
enough to make them observable. The experimental technique most widely used to study trions is time resolved 
photoluminescencc 14,15 . Processes as generation and recombination of excitons and trions have been recently analyzed 
through this technique^. 
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FIG. 1: Scheme of the asymmetric double quantum well structure under free electron resonant condition. 
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In this work we propose a new method to generate trions in coupled DQW through the direct-created excitons and 
the necessary free electron excess. In the usual experimental conditions, donor concentration is greater than acceptor 
concentration (~ 10 18 and ~ 10 16 cm -3 , respectively). The process of trion generation can be considered in three 
steps: first, if the Fermi level of the n-type material (see Fig. 1) resonates with the electronic level in the left-hand 
wide quantum well (LQW), electrons tunnel to LQW forming the cross or spatially indirect (interwell) excitons with 
holes diffused into the narrow right-hand well (RQW). Second, if the resonance condition for electrons in the excitonic 
state between both wells is satisfied, electrons in the excitonic state tunnel to the RQW forming spatially direct 
(intrawell) excitons. Third, if there is an excess of electrons with respect to holes, free electrons will coexist together 
with direct-created excitons. Some of these free electrons can be bound to excitons leading to trions. It should be 
noted that, in coupled quantum wells, the strong dipole-dipole repulsion between excitons prevents the formation of 
biexcitons. Thus, we will neglect the contributions from these neutral species^. 

Under resonance conditions the charge density in coupled DQWs oscillates between both wells leading to an os- 
cillating dipole moment and the corresponding Terahertz radiation emission. If we consider the holes confined in 
the well where they were initially injected (due to their large effective masses), the electronic charge will play the 
main role in the coherent process. Because we will deal with three different types of electron binding (free electron, 
excitonic electron and trionic electron) they will occupy different energy levels. Moreover, the resonance conditions 
will be slightly different for each case and thus, their oscillatory behavior will be different as well. 

The aim of this paper is the theoretical analysis of the electron in trion dynamics in presence of direct-created 
excitons and free electrons. We will carry out this in the framework of the matrix density formalism, taking into 
account the different generation, recombination and annihilation rates for electrons with different type of binding. 
The dynamics of coexisting excitons and free electrons was studied elsewhere^. The work is organized as follows: in 
Sec. I we present the exciton and trion Hamiltonians and the corresponding Schrodinger equations, as well as the 
dynamics of the electronic density evolution. In Sec. II we include numerical results for the charge density dynamics 
and for the current density. We also include some comments regarding the results and the used approximations. 

II. FORMALISM 

Initially, a group of electrons tunnels into the LQW. The charge, dynamically trapped within the wells, produces a 
reaction field that modifies the structure profile and then, the time evolution of the system. Thus, we need to calculate 
wave functions for the free-electron, the excitonic electron and the trionic electron, and the corresponding electronic 
levels, by solving self-consistently the one-particle Schrodinger equation, together with the Poisson equation for each 
species. 

Several techniques have been used to solve the Schrodinger equation for neutral excitons and charged excitons 
(trions). We can cite the commutation technique, developed by Combescot et al^. Although it is a very attractive 
formalism, it is almost inapplicable from the numerical point of view. Another technique, based on the deformed 
correlated Gaussian wave functions, is the Stochastic variational method used by Riva et al^S. Whittaker and Shields 20 
used a set of Landau-level basis to study the in-plane (xy) motion. These two methods are specially indicated when 
an external magnetic field is applied. In presence of longitudinal electric fields and in absence of magnetic fields 
Dacal and Brum2i have developed a technique based on the configuration interaction method ( CIM) 22 . Because of 
its feasibility we will use this last technique with some implementations. 

The CIM considers that the longitudinal (z) and transversal (xy) parts of the wave functions are separable. Thus, 
we first factorize wave functions of excitonic and trionic electrons into the corresponding in-plane part and confinement 
(growth) direction part. We can't neglect the in-plane transversal motion part of the Hamiltonian because both parts 
of the split Hamiltonian are connected by the Coulomb potential. Such a circumstance affects wave functions of the 
different types of electrons. In order to reach the resonance for the different electronic states we have included in 
calculations an external electric field applied in the growth direction. The confinement part of the Hamiltonian can be 
solved by means of the Airy functions within the effective mass approach for electrons and holes. To solve the in-plane 
part of the one particle Hamiltonian we will follow the same steps as in Refi^i. First we will use polar coordinates in 
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terms of the center of mass (CM) and relative to the hole coordinates. To say, 
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Phi 
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m h ±P h 



m h± + 2m e 

where nih± is the in-plane hole effective mass. We have considered an isotropic electron effective mass. 



(1) 



A. Exciton and trion Hamiltonians 

The excitonic Hamiltonian can be written as 

H exc = H(z e ) + H(z h ) + H ± + V c , 

where 

h 2 d 1 d 



H(z e ,h) 



2dz e h m e h z dz e h 



V D Qw{z e .h) + V H {Ze,h) ± |e| Fz e . h . 



(2) 



(3) 



Here m^z is the hole effective mass in the growth direction, VDQw{ z e,h) is the asymmetric DCW potential, F is the 
electric field required to reach the inter-level resonance. VH{z e ,h) is the Hartree potential caused by the Coulomb 
interaction of the electrons with themselves and with the donors (or holes with acceptors). The self-consistent potential 
VH{z e ,h) can be obtained by solving the Poisson's equation 



d 2 V H {z e . h ) _ 47re 2 



N D ,A{z e ,h) ~ \<l>j( z e,h)\ 



(4) 



where Nd,a is the doping profile; thus, / ND,A(z e )dz e ,h = n je,h = n 2 De,h is the in-plane averaged 2D e ^ h density 



of electrons (donors) or holes (acceptors) and j refers to the occupied subbands. 
The in-plane term reads 



2/i 



1 d 



d 



p dp \ dp J p 2 dp 2 



i d 2 



(5) 



where we have considered invariant the in-plane excitonic p reduced mass. The Coulomb potential is 



Vce,h(P) 



dz e dz elt h4>l(Ze)4>l'h( z e',h)- 



P 2 + (Z e - Ze',hY 



(6) 



and includes electron-hole (e — h) and electron-electron (e — e') interactions. This Coulomb effective potential is valid 
for excitons and for trions, being slightly different for each case due to their different electron wave functions and 
expected values of (z — z'). We will assume a constant dielectric function e across the interfaces. For practical reasons 
we change the notation of the above Coulomb potential to 



v ij ( P ) = e -^ 



dzdz' 



yJfP + (z-zf' 



(7) 



where ej = — e for i = u,l and ej = e (— e) for j = h (u,l). Indexes u, l,h refer to electrons in the upper or lower 
levels and to holes in the right-hand QW, respectively. 



4 



The trionic Hamiltonian can be written in an analogous way to the excitonic case by using relative coordinates for 
the transversal motion: 



H tr = V [H(z ei )+H X i]+H(z h ) + ^— pip 



1=1,2 
^2 



m h ± 



Pi ~ P2I + ( Z el ~ Z e2 ) 2 



(8) 



where is the in-plane linear momentum for the ith particle, 

h 2 d 1 d 



H{z ei . h ) = 



2 dz e i : h Tti e iji Z dz e i h 



V D Qw(z eith ) ± |e| Fz ei , h , 



and 



2/i 



1A f A^i 1 92 

Pi dp, 1 ' J pf dp 2 



eJ (z el - z h ) + p 2 



(9) 



(10) 



To solve the above Hamiltonians we will follow the same steps as RefiSi. Thus, through the CIM we construct a 
non-orthogonal basis set of Slater determinants to solve the eigenvalue problem. This method is widely described in 
Refi22. We build up the exciton trial wave function as 



ijk 



(11) 



where Cijk is a variational parameter, Nijk is the determinant normalization, indexes ijk refers to one particle i, j, 
k, respectively, p — p e — p h , m is the relative particle angular momentum, and 4>(z e ^h) is the envelope wave function, 
and </9™(p) is the one particle in-plane wave function 



ipf(p) =N kmP m exp 



A 



exp (imO) 



(12) 



Here Nk m is the relative particle normalization factor and X k is a set of parameters (one for each angular momentum 
m) that we have taken from RefiS. 

The above expressions can be generalized for trions as 



ty tr (z e l, Z e 2, Zh, Pl, P 2 ) 



E ■ 

ijmnpqr 



'ijmnpqr -^ijmnpq 



4> q {Zel)tpf (Pi) 4> r {Ze2)^ (p 2 ) 
+ r (z e i V™ ( Pl ) <j) q {z e2 )ipT (p 2 ) 



(13) 



It is important to remark that, in absence of external magnetic fields, the only bound state for trions is the singlet 
to + n = 0. We would like to mention here that we have included a new contribution: the spatial bending of the 
structure due to the doping. Thus, after variationally obtaining the above wave functions for VH{z e ^h) = 0, we need 
to solve the Poisson equations (4) and turn again to the Slater Hamiltonians. Fortunately, we only need to modify one 
particle envelope functions in the growth direction by solving self-consistently the z-direction part of the corresponding 
Schrodinger equations (3, 9) together with the Poisson equations. For practical purposes, to numerically work the 
z-direction part with the Transfer Matrix Method is better than using the Airy functions. Variational coefficients and 
normalization functions for wave functions (11, 13) will be slightly different when including the spacial self-consistent 
potential Vn{z e ^). Based on Ref^ we use parameters for the in-plane part of the Hamiltonians corresponding to 
s and p ± states. The inclusion of the cF states only affects about a 5% the position of the excitonic and trionic energy 
levels but it supposes an additional computational effort which is not reflected on the final results. Another key point 
consists on the identification of the levels because it will be fundamental in the next section. This information is 
contained in the envelope functions (f>(z e ^h)- Thus, both for electrons and holes, the two deepest levels correspond to 
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FIG. 2: Energy levels versus applied electric fiel. 
trionic electrons. 
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Dashed line: excitonic electrons. Dotted line: 



the resonant levels of the two coupled wells. Hereinafter we will distinguish them with the indexes u (upper) and I 
(lower). We have neglected the possible hole states mixing for the applied electric fields under consideration because 
we are mainly interested in the conduction band states. 

Concerning to free-electron, the electronic levels (29.97 meV and 36.61 meV) would correspond to the LQW and 
the RQW, respectively, if they were decoupled. For the excitonic electron case these two levels are shifted by the 
Coulomb interaction by different energy amounts corresponding to the indirect and direct excitonifl ( 5.78 meV and 
8.60 meV, respectively). We have considered that hole energy levels remain the same due to the bigger hole effective 
mass. The case of trions is more complicated due to two possible additional forms of setting the new electron with 
regard to the previous excitonic state. From the energy point of view, the most favorable cases correspond to the 
direct exciton bounded to an electron either in the LQW or in the RQW. That means a shift of the free electron 
levels of 6.91 meV and 9.93 meV, respectively. We represent in Fig. 2 electronic levels as function of the external 
electric field. It can be seen that the resonances of the three species take place for almost the same applied electric 
field. With these values we calculate the level splitting for coupled wells, Ax, and the transmiss ion matrix element, 
T, which equals to Ay/2 at resonance. The level splitting energy for decoupled wells is A = \/A§, — 4T 2 , for each 
case. We will use these parameters in the next subsection. 



B. Time evolution of the electron density 



In this Section we will treat the dynamic behavior of the electronic density within the Matrix Density formalism 
in the momentum representation. Although electrons spend some time in tunneling the left-hand confining barrier, 
we assume an ultrafast S(t) injection of electrons in the LQW to simplify calculations. Also, diffusion of holes from 
the p-doped material to the RQW is considered as a 8(t) function. We use this approximation because we are mainly 
interested in the dynamics of the electronic density between wells after injection. Based on the previous premises, we 
start from the general quantum kinetic equation for the density operator R(t) 



dR(t) 
dt 



H, R(t) 



0. 



(14) 
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Here H is the DQW many-band Hamiltonian we have described above. If we project on the conduction band states 
we find the 2x2 matrix kinetic equation in the momentum representation 



d [p P (t)} 



or 



at 

G P (t) 



^C' a j'sp a O s P 



00' 



(15) 



This equation describes the time evolution of the density matrix p p (i), where indexes j, j' refer to the levels u (upper) 
and I (lower) of the DQW, and p is the electronic wave vector. In the expression (15) [/3 p (£)]_.,., = ^at'sp a J s p)> 
where cij, and a,j sp are the creation and annihilation electron operators, s is the electronic spin, and ( ) means the 

statistical average. The operator h = (A/2)a z + Ta x is the matrix of the one-electron DQW Hamiltonian without 
Coulomb interactions, and 3^ (i = x, y, z) are the standard Pauli matrices. For numerical and analytical reasons it is 
convenient to remove the Coulomb interaction from Hamiltonians (3) and (10) and to work it separately in the Hartree 
representation. Thus, variational coefficients and normalization factors in 4>™(p) (13) will be slightly different. The 
Coulomb interaction can be written as 



isk+q"js'k'-q"JS 



(16) 



ij ss' kk'q 



where i, j refer both to electron and hole states to include electron-electron, electron-hole and hole- hole Coulomb 
interactions. S is the normalization area and q = |q|. The matrix elements Mij(q) are defined as^ 



M l3 {q) 



2neiej 
e<7 



dzdz' exp (—q \z — z'\) ^(z)! 2 U-(z')| 



(17) 



where — — e for electrons = e for holes (e is the electron charge). 
After substitution of (16) in (15) we reach 

dp P (t) i_ 

at n i 



h + Q(t)+X p (t),p p (t) =G p (t), 



(18) 



where Q(t) is the direct Coulomb potential and X p (t) is the exchange potential. If we define the projectors over the 
I and u conduction states as Pi = (1 + <r z ) and P u = (1 — a z ), respectively, Uh{t) is the hole concentration in the 
h-th hole subband, and ni. u (t) — (2/ S)J2 p Tr Pi tU (tJp(t) is the electron concentration in the upper (u) and lower 
(I) levels, thus matrices Q(t) and X p (t) can be written as 



Q(t) 



Mu(q)m(t) + Mi u {q)n u (t) ~^M lh {q)n h {t) 



P 



M uu {q)n u {t) + M u i(q)m(t) - ^ Mi h (q)n h {t) 



Pu 



(19) 



| , [M u {q) - Mi u {q)]Tr [Pip^t)) Pi 

X p^ = s E \ + [M uu (q) - M lu {q)] Tr (i 5 tl p p _ q (t) 

+M lu (q)p p _ tl {t) 



(20) 



After integrating (18) over momentum p one finds that the non-diagonal part of the exchange potential , 
Mi u (qyp p _„(t), does not contributes to the commutator in eq. (18). The diagonal part of X p (t), together with 
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the Coulomb potential Q(t), provides us the renormalization of the levels and, thus, the renormalization of the split- 
ting A because these local terms couple different subbands. It is important to point out that, contrary to Ref. (3), 
in our case we are dealing with typical values of q of the order of 10 s m _1 , which is the same order as the distance 
between the centers of the wells Z. Thus, we cannot neglect the exchange contribution. 

Now, we define the concentration matrix as n(t) = (2/5) Pp (*) an d integrate over p to obtain the isospin 
representation^. - 



n(t)=n (t)+J2n i (t)a i {t)=no(t)+n{t)a. 



(21) 



After substituting in (18) and introducing the subindex k to differ free electron (/), excitonic electron (exc) and 
trionic electron (tr), the above Liouville (18) equation leads to the Bloch-equation system for the evolution of the 
isospin density vector n./-(t) = (n%(i) , rit(t) , ?if(i)) and the scalar density n k (t) : 



dn k (t) 
dt 



- [L k {t) x n fc (t)] = G fc (t)-A fc (t), 
9n°(f) 



dt 



N k S{t) - S k (t). 



(22) 



where N k 6(t) - S k (t) = (l/S)J2 p Tr (d pk (t)) and G k (t) - A k {t) = (l/S)J2Tr(aG pk (t). The nonlinear dynamic 

properties of the system are described by L k (t) = (2T k /h, 0, A k (t)/h), where A k (t) is the renormalized time dependent 
level splitting induced by the carrier density and T k is the tunneling matrix element. G k (t) = (0,0, N k S(t)) is the 
generation term and N k denotes the initial density of electrons. In the present case, just at t = 0, Nt = N e (injected 
electron density) and N exc , tr = 0. The terms A k (t) and S k (t) include generation and annihilation rates for excitons 
and trions. We will describe them below. Vector equations (22) split into three coupled Bloch system, corresponding 
to free electrons, excitonic electrons and trionic electrons. The time dependent energy splitting renormalization is 
determined (after an easy but tedious algebra) by the many-body Coulomb interaction through the Hartre-Fock 
approximation (19,20) 



A exc (tr)(t) — A e 



2ire 2 



a(tr) 



dzdz \z — 



u,exc(tr) 



(*') 



■^l,exc(tr) 



u,exc(tr) 



CO 



i u u ^i) - 



l,exc(tr) 



c(tr) 



(t) 



l h (t) 



27re 2 



dzdz'8(z — z') 



l,exc(tr 



.(*') 



u,exc(tr) 



(z') 



H,exc(tr)\Z) 
^M,exc(tr)( Z ) 



exc(tr) 



(0 



c(tr) 



(23) 
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and 



A f (t) = A f + 



27re 2 



dzdz' \ z — z'\ 



K/C^r-i^,/^) 

\<t> uJ (z)\\}(t) + \^ u (z)\ 2 n l f (t) 
2ne 2 



dzdz'S(z-z') \(f> lif (z')\ - \(j> u j(z')\ 



(24) 



where (i) and fi/i(i) are the electronic density in the upper (lower) level and the hole density in the right-hand 
well, respectively. The last terms in (23) and (24), proportional to the coefficient a, correspond to the exchange inter- 
action contribution. To calculate this contribution we need to introduce the carrier screening since the use of the bare 
exchange potential may lead to divergent results due to its singularity when q — > 0. This means that the energy renor- 
malization caused by the exchange term may diverge. In an easy approach, and considering parabolic subbands filled by 
degenerate electrons, we use the statically screened Coulomb potential V s (q) = e (^ q ) i where q s = ^p- is the screen- 
ing wavevector. In this case, the screening exchange can be written as y)V^(q) [p(p — q)]„„ +4lim [V^M — ^( r )]- 

The last term corresponds to the difference between the screened and unscreened Coulomb potential, in the real 
space, at r = 0. Integration over p and q of the exchange commutator in (18) leads to corrections of the order of 
11% in the level splitting. In the present work a « 12. 1A. Exact calculations lead to three different values of the 
parameter a, but differences (of about 10%) are not visible in results. 

We will assume that only the first heavy hole level is occupied. It is worth mentioning that excitonic and trionic 
electrons not only interact among themselves and with holes but with free electrons as well. It is important to 
remember that, for direct-created excitons in DQWs, holes and electrons are generated in different wells. Equations 
(23, 24) differ from similar expression in Ref^ due to the exchange contribution. 

If we denote by Nh as the initial total hole density, then nh (i = 0) = Nh and 



n(t) = 4(t) + "L(t) 



n&(t) + n\ r {t) 



(25) 



. The factor 1/2 in the last addendum appears because n\ r denotes the trionic electron density which is twice the 
hole density forming trions. The scalar density n°(i) is the sheet electron density (the total density of electrons in 
the DQW structure). Also, we can write 



nl(t) 



n l k {t)-nt{t). 



(26) 



Now we define r as the relaxation time for the non-diagonal part of the density matrbsiS^. This relaxation leads to 
the interlevel redistribution of the electron density. We consider here the generation and annihilation terms. Initially, 
free electrons are injected quasi-instantaneously and lose their free condition to generate excitons and trions. So we 
have two annihilation mechanisms expressed through the generation rates T exc and r tr , respectively. Moreover, trions 
destruction gives rise to free electrons because the recombination of one electron with the hole leaves an additional 
free electron. Thus, we have an additional generation term for free electrons at a rate l/rt r . By writing densities 
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rif(t) and n l exc tr (t) in units of N e , equation (22) for free electrons becomes 

l„»( t )-^„5(t) + ^n3l(t) + inj(*) = 0, 
-n}(7i) - -^n y f (t) - -fS{t) n° tr (t) + T tr n m (t) + T exc n h (t) = 0, 



dt 1 W ft J W N e ^ ' T tr 

^n%t)-^5{t)-—n%(t)+T tr n m {t)+T exc n h {t) = 0, (27) 

at ' N P Ttr 



where n m (t) = min[nj(t),n® xc (t)]. In a similar way excitonic electrons appear, after the instantaneous hole diffusion, 
at a rate T exc and disappear due to trion generation and recombination at rates T tr and l/r exc , respectively. Thus, 
for excitonic electrons, 

J f < xc (t) + ^ftnl xc {t) + In^t) = 0, 

jn y exc {t) - ^Y^< xc {t) + ^pnl xc (t) + \nl xc (t) = 0, 

A 2T 1 

(t) - ~^nl xc {t) - T exc n h (t) + n° exc (t) + T tr n m (t) = 0, 



71 

dt excy "' h 



and, for trionic electrons, 



^n° exc {t)-T exc n h (t) + — n a exc (t)+T tr n m (t) = 0, (28) 



|nf r («) + 





1 X 


(*) 


= o, 


+ <(t) 


1 V 


(*) 


- o, 


- r tr n m (t) + 


n tr 

T tr 


(*) 


= o, 


- r tr n m (t) + 


n tr 

T tr 


(*) 


= 0. 



d z /.\ ^Ttr y 

-n tr (t) - —nl 

±n° tr {t) - T tr n m (t) + —nl(t) = 0. (29) 

dt Ttr 

In the above expressions T^ 1 is the mean time required for the trion formation and T exc and r tr are the intrinsic 
exciton and trion lifetimes. An extensive analysis about these times can be found in the paper of Esser et al^. We 
have considered that the process of exciton formation is much faster than the generation of trions and the annihilation 
of both excitons and trions^. Actually, T tr corresponds to the scattering rate of excitons with free electrons, which 
gives rise to the trion formation. Experimentally, this scattering formation rate is proportional to the free electron 
density, to the in-plane exciton area (which contributes to the cross section of the process) , and to the kinetic energy 
of the two interacting species. In our work we have used data for GaAs — GaAlAs included ini& . 
Equations (28), (29), and (30) allow us to get the total density values 

n z T {t) = n z f (t) + n z exc {t) + n z tr (t), 

n°(f) = n f (t) + n° exc (t)+n° tr (t), (30) 
where n^,(t) is proportional to the total dipole moment. 

C. The current density 

Within the dipole approximation the current density can be written as 

3w(t) = jE^(*>i«" ( 31 ) 
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where I, I' refer to the coupled levels, v±_iv — {4>i( z )\ -§z \4>i'{ z )) — v J-iu = — ^ {<l>i(z)\ -§z I0 U ( Z )) is the intersubband 
transverse velocity, and L is the size of the structure. Since we have supposed a left-hand barrier doping to locate the 
Fermi level just above the two deepest conduction levels, only these two levels will be occupied and we can neglect 
contributions from remote bands. Therefore, the whole problem is now reduced to the additional calculation of the 
transverse velocity. 

III. RESULTS AND DISCUSSION 

We have used in calculations an asymmetric DQW (Fig.l) configured by a llOA GaAs LQW and a 100A GaAs 
RQW, separated by a 30A Gao.55Alo.35As barrier. Then we can guarantee the interaction of the electron wave 
functions of neighboring wells (non-null overlap). For this structure the distance between the centers of the wells is 
Z c = 135A. We have also considered the permittivity as constant along the whole structure, e = 12.9. To analyze 
the contribution of trions to the total dipole moment we perform numerically the coupled system (28-30) through 
the Runge-Kutta method. We will restrict ourselves to cases in which the density of injected electrons is higher than 
the density of diffused holes, e.g., an excess of free electron density of about 10 9 — 10 10 cm~ 2 . Higher densities could 
prevent the trion formation^. 

Figures 3(a-c) represent the density evolution of free, excitonic and trionic electrons [«.°(i)] for different Nh/N e 
ratios. Initially, there are only free electrons. The density of these electrons decays quickly, giving rise to the 
excitons. These excitons yield to trions generation due to the scattering processes between excitons and free electrons. 
The recombination of electron-hole pairs reduces the density of electrons both in exciton and in trion. In the last 
case each extinct trion leaves one electron free, which contributes to increase the stock of remaining free electrons. 
That recombination leads to a long time limit of the total density of electrons that equals the free electron excess, 
n^(t — > 00) = n°j(t — » 00) = 1 — Nh/N e . Cases (a), (b), (c) corresponds to Nh/N e — 0.8, 0.5, 0.2, respectively. At 
first sight, the formation of trions seems to be more efficient when the initial hole concentration is half the initial free 
electron concentration. In this case the trionic electron density can even exceed the excitonic electron density. 

Next we analyze the behavior of nj*(i) to obtain n^(t), which is proportional to the dipole moment. As stated in 
Eq. (26) n^{t) is the difference between the densities of the lower and the upper energy levels. We represent n z k (£) 
in Panels 4(a — c) for N e = 10 10 cm~ 2 . Each figure corresponds to one of the above mentioned Nh/N e ratios. As 
expected for low densities, free electron density oscillates between lower and upper levels with an exponential-like 
envelope profile not only due to the depopulation of the corresponding levels but also to the interlevel redistribution. 
However, Coulomb interaction forces the excitonic electrons to oscillate staying basically in the upper level before 
the recombination and the interlevel redistribution. Trionic electron tends to occupy mainly the lower level. On the 
whole, the mentioned contributions give place to a modulation of the free electron dipole moment which diminishes 
when Nh/N e does. 

For higher densities the effect of the electron-electron Coulomb interaction becomes more noticeable. Panels 5(a — c) 
represent n^(t) for N e = 5 x 10 10 cm~ 2 . Free electron density behavior is absolutely different than the low density 
case and, initially, most of the electrons stay in the lower level and only a few of them oscillates between wells. At 
the same time excitonic electrons experience a strong confinement in the upper level and trionic electrons do the 
same in the lower level. The frequency of the three species of electrons shows an increase that is more remarkable 
for the trion case. As in the low density case, the effect over the free electron dipole moment due to excitons and 
trions decreases when Nh/N e decreases. Another effect of the charge density increase is to diminish the oscillation 
period breaking down the harmonic response of the dipole moment. The anharmonic behavior can be explained as 
a consequence of the space-charge potential created by the spatial distribution of electrons and holes. This space- 
charge potential is repulsive for holes and attractive for electrons. For densities higher than the values used in this 
work (N e > 10 11 cm~ 2 ) the Hartree-Fock potential becomes more important and a better approximation to analyze 
the carrier distribution would be necessary. However, the Hartree-Fock approximation used here for the many-body 
Coulomb potential incorporates the main peculiarities of the particle interactions as function of the carrier densities. 
Thus, electron-hole attraction dominates over electron-electron and hole-hole repulsion at low densities. When the 
carrier density increases the repulsive part of the Hartree-Fock potential increases as well. Beyond certain initial 
density of free electrons, the repulsion equals the attractive Coulomb interaction and the oscillation period decreases 
for the three species. For higher densities the Coulomb interaction of the second electron with the hole in the trion 
is canceled resulting in the trion extinction, which leads to excitons and free electrons. If the free electron density 
increases even more, the binding energy of excitons tends to zero and excitons also disappear. 
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FIG. 3: Temporal evolution of the normalized electron density, (a) N h /N e = 0.8, (b) N h /N e = 0.5, and (c) N h /N e = 0.2. Solid 
line: free electrons, dashed line: excitons and dotted line: trions. 



The shape of the current density is the same as the shape of the n^(t) density and we don't believe necessary to 
reproduce them here. As an example we include current quantum beats in Fig. 6(a) where upper panel corresponds 
to the current density for parameters of Fig. 4(a), and lower panel to data of Fig. 5(a). It should be notice that, in 
any case, current density vanishes for long times due to relaxation as Fig. 6(6) shows for the same parameters as in 
Fig. 6(a). 

Now, we will mention here the approximations used in the present method of calculation. We assumed that the 
electron injection and hole diffusion are processes much shorter than the coherent oscillation period and the exciton 
generation time, which, in turn, is even shorter than trion generation and recombination times. Another important 
point is the possible effect of the applied external electric field on the trions. We have used field intensities lower than 
5kV/cm to avoid the ionization or the possible trion diffusion along the structure^. 

In summary, we propose a new method for the generation of trions that avoids the undesirable effects coming from 
the interaction of the laser electromagnetic field with photoexcited carriers, an ever-present problem in photoexcitation. 
The method is based on electron tunneling injection and hole diffusion from remotely doped layers^. We have also 
analyzed the effect of these trions, together with the direct-created excitons, on the time dependent dipole moment 
caused by the free electron density oscillations. Such an effect is only relevant for an initial excess of the free electron 
density with regard to the hole density. We have studied the above mentioned effect through the three coupled Bloch 
system obtained from the Liouville equation and taking into account five possible time dependent processes. Our 
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results show a wide variety of responses, from modulation and quantum beats to anharmonic regime, caused by the 
different carrier densities of the three types of electrons. These densities mainly affect the dynamics of the level 
splitting energy of electrons and thus, the coupling strength. The combined electron oscillations yield to a nonlinear 
coherent electromagnetic radiation from the semiconductor structure. In other words, the existence of trions could 
be detected by means of the current density, which is straightforwardly translated into Terahertz radiation. 

As far as we know only photoluminesccncc experiments for trions in coupled DQW are available. Moreover, there 
are not any studies, neither experimental nor theoretical, about the possible effect of trions on the dipolar radiation 
emission from quantum wells. Therefore, we are forced to extrapolate from well-known results for excitons. Since 
the standard way to generate trions is ultra-fast photoexcitation, a detailed analysis of this process deserves special 
attention together with the peculiarities of the high-density case. We expect that the presented work will stimulate 
experimental studies of the dynamics of direct-created excitons and trions. 
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FIG. 5: Temporal evolution of the normalized dipole moment with N e = 5 x 10 10 cm~ 2 . Upper panels: free and excitonic 
electrons. Lower panels: trionic electrons and total dipole moment, (a) Nh/N e = 0.8, (b) Nh/N e = 0.5, and (c) Nh/N e = 0.2. 
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